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We propose and support the possibility that the shape of topological density 2-point function in 
pure-glue QCD is crucially, and possibly entirely, determined by the space-time folding (geome- 
try) of the double-sheet sign-coherent structure of Ref. [1], while the distribution of topological 
density within individual sheets only determines the overall magnitude of the correlator at finite 
physical distances. A specific manifestation of this, discussed here, is that the shape of the cor- 
relation function (encoding e.g. the masses of pseudoscalar glueballs) is reproduced upon the 
replacement q(x) — » sgn(q(x)), i.e. by considering the double sheet of the same space-time ge- 
ometry but with constant magnitude of topological density. Combined with previous results on 
the fundamental topological structure, this suggests that a collective degree of freedom describ- 
ing topological fluctuations of QCD vacuum can be viewed as a global space-filling homogeneous 
double membrane. Selected possibilities for practical uses of this are discussed. 
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1. The Context. Discovery of sign-coherent topological structure in equilibrium configurations 
of lattice-regularized QCD [[Tj] and the subsequent demonstration that this order is of dynamical 
origin opened the door for systematic inquiry into the nature of configurations dominating the 
QCD path integral [J|]. Conceptual innovations associated with these developments have their 
origin in the fact that the space-time structure is detected directly in equilibrium ensembles and 
using local composite operators which, among other things, has two important implications. (1) 
The information on the space-time order in QCD vacuum so obtained is free of a priori assumptions 
and subjective manipulations {Bottom-Up approach [||]) thus putting the associated line of research 
on a solid ground. (2) The space-time structure identified in this way incorporates features at all 
scales in the continuum limit. Consequently, the association of these features with underlying 
physics is not limited to usual "vacuum structure" low energy manifestations (condensates, string 
tension) but extends to arbitrary property described by the theory (fundamental structure) ^ ^J. 

In the general framework of ^ that we follow, it is emphasized that the scale dependence of 
phenomena in QCD has to be reflected in the way we describe the structure of typical configu- 
rations (vacuum structure). Indeed, while the fundamental structure carries complete information 
on the theory, it is implicitly understood that we are simultaneously considering an infinite tower 
of effective structures, labeled by momentum scale A. In the effective structure the fluctuations 
up to length scale 1/A are averaged out from the fundamental structure [|3j]. An individual ef- 
fective structure does not carry complete information about the theory, but the association of its 
space-time features with physics at the corresponding scale A is expected to capture this physics 
most efficiently. Loosely speaking, effective structures represent the fundamental order at varying 
resolutions (scales of physics). Here we focus on the fundamental structure in topological density. 

The chief rationale of inquiries into the structure of typical QCD configurations is the expec- 
tation that such information can eventually be turned into an improved qualitative and quantitative 
control over the theory. In case of fundamental topological structure, the underlying order man- 
ifests itself via existence of extended sign-coherent lower-dimensional regions ("sheets") with 
strong spatial correlation between the regions of opposite sign [jj]]. More precisely, the follow- 
ing properties of the sign-coherent structure have been advocated, (i) Lower dimensionality [fil], 
meaning that it is impossible (in average sense) to embed a 4-d sign-coherent ball of finite physical 
radius into the sign-coherent topological structure of QCD (see also 0, |5|]). (ii) Inherent glob- 
ality [JXJ, [(J, meaning that it is not possible to break the coherent structure into localized pieces 
without compromising the physics. In fact, there are only two global oppositely charged sheets 
closely following each other ("double-sheet"). (Hi) Space-filling nature meaning that even 
though the structure is locally lower-dimensional, it nevertheless fills a macroscopic fraction of 4-d 
space-time and is thus geometrically analogous to space-filling curves. (For additional references 
related to fundamental topological structure, see [{7|].) 

In this work, we propose a new geometric property of the double-sheet structure that has direct 
connection to QCD dynamics. In particular, we will show that the shape of the 2-point function 
is insensitive to inhomogeneities within the sign-coherent regions present in typical configurations 
of lattice-regularized ensembles. In other words, the presence of such inhomogeneities is due to 
an unphysical noise, not affecting correlations at finite physical distances. Rather, at the level of 
fundamental structure, the dynamics appears to be encoded (possibly entirely) in the space-time 
folding of homogeneous objects - in the geometry of the double-sheet sign-coherent structure. 
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Figure 1: The results for the original, G(r), and the geometric, G(r), correlators on 24 4 lattice in pure-glue 
lattice gauge theory (Iwasaki action) at a = 0.055 fm. G(r) was rescaled to coincide with G(r) at r = 3. 

2. The Basic Observation. Original results put forward in this work are based on the outcome of 
the following numerical experiment. Consider an ensemble of topological density configurations 
corresponding to Wilson's lattice gauge theory at some cutoff. With every configuration ^ = 
{q(x), x G Z 4 } associate a new configuration ) defined via 

^={ sq ( x )= sg a(q(x)),q(x)e^} (1) 

where sgn(y) is the sign function. Thus, the configuration ^ has field values ±1 and possesses 
the same sign-coherent regions as c €. Now, consider the two-point functions of the original 
G(x) = (q(x)q(0)) and of the associated G(x) = (sq(x)sq(0)} ensembles. Is there a definite re- 
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lation between G{x) and G[x)l Note that we will frequently refer to barred entities as "geometric" 
(i.e. geometric correlator, geometric ensemble) since the only information kept from the original 
configuration is the space-time geometry of sign-coherent regions. 

In Fig. we show the result of such calculation on 24 4 lattice at lattice spacing a = 0.055 
fm, using overlap-based definition of topological density. Since the overall magnitudes of the two 
correlators will evidently be very different, we have rescaled the geometric one such that at \x\ = 
r = 3 they become equal, i.e. G(r) — > G(r)G(3)/G(3). On the bottom plot, we show the results 
on a logarithmic scale to better discern points at larger distances. As can be seen quite clearly, the 
original and the rescaled geometric 2-point functions appear identical within statistical errors. In 
fact the agreement is better than statistical, indicating that, already at this physical volume, sign- 
coherent geometry drives correlations at the level of individual configurations. 

There are couple of points that we would like to note here. Firstly, the feature uncovered above 
only holds in the "negative part" of the lattice correlator, i.e. in the region where it can approximate 
the behavior at finite physical distances, and not in the positive core. Indeed, we will not propose 
that lattice correlators G(x) and G{x) are equal up to rescaling. Rather, we will conjecture the 
possibility that they define the same continuum limit at non-zero physical distances. Secondly, 
we emphasize that the currently available data [^] is supportive of the scenario with geometric 
dominance being satisfied at all, rather than just at asymptotically large physical distances. 
3. Specifics and Formalization. The numerical work on which we draw our conclusions is based 
on pure-glue SU(3) lattice ensembles defined by the Iwasaki action and specified in Table |l| The 
referenced lattice spacings are obtained from string tension (yfc =450 MeV), and the physical 
volume is kept fixed at V p = 3 fm 4 . For definition of topological density operator we use overlap 
Dirac matrix [§] based on Wilson-Dirac kernel with standardly defined parameters r = 1 and p = 
26/19. This setup was used in the original work [Q] as well as in all the followups by the Kentucky 
group. Implementation of matrix-vector operation needed to evaluate topological density q(x) = 
j^try5D(x,x) is described in Ref. [|l0|]. 

The correlators shown in Fig. [I] are coarse-grained over the distance of half a lattice spacing. 
By this we mean that the samples for correlations with a given point were accumulated within 
the spherical shells of thickness one half, centered at that point. More precisely, let us define a 
sequence of points and the associated sets (intervals), namely 

_ f for* = _ f {0} for* = 

n \ k 2 + \ fbr* = l,2,... Jk \[r k -lr k + \) for* = l,2,... {) 

The coarse-grained correlation function Gfa,^) on a given configuration ^ is defined as 

G(r k ^) = (q(x)q(y))^ ]x _ y]eJk (3) 

Continuum limits of the (ensemble-averaged) standard and coarse-grained correlators are the 
same, but the latter behaves more smoothly which is why we use it here. All conclusions pre- 
sented here are independent of this choice. In what follows we will denote the coarse-grained 
correlator simply as G{r) with the above discrete values of r implicitly understood. 

The results shown in Fig. [j] correspond to the ensemble $4 with the ensembles at coarser lattice 
spacings behaving similarly. To quantify the trends and to formalize the proposed observation, we 
need a suitable measure. Equality of functions up to a multiplicative constant is naturally expressed 
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ensemble 


a [fm] 


V 


V p [fm 4 ] 


configs 


£\ 


0.110 


12 4 


3.0 


142 




0.082 


16 4 


3.0 


50 


^3 


0.066 


20 4 


3.0 


20 


£4 


0.055 


24 4 


3.0 


10 



Table 1: Ensembles of Iwasaki gauge configurations for topological density calculations. 

in terms of a normalized overlap (scalar product). For continuum real-valued functions f\ , f2 on 
the interval [a, b] one has 

0\f h f 2 ,a,b] = fi-f 2 = I dxMx)f 2 (x) (4) 

|/l|l/2| J a 

but in case of fast-decaying functions (such as the correlators considered here) a measure that is 
more uniform over the domain and more stringent is provided by the "relative overlap", namely 

0\fi,f 2 ,a,b] = 0[l,^,a,b) (5) 

Note that for functions f(x k ) defined only on discrete set of arguments x k (lattice correlators) we 
apply the above formulas as well via replacing these functions with their piecewise constant exten- 
sions f{x) = f{x k ) for x e [x k ,x k+ i). 

Consider the lattice definition of the topological density 2-point function in the domain of 
physical distances [r p ^,r p ^\- This requires studying the behavior of G(r,a) on the ever-growing 
sliding lattice interval [r p j/c?, r p ^/a] close to the continuum limit a — ► 0. If the original and the 
geometric correlators (G(r,a), G{r,a)) carry the same information about the shape of the physical 
2-point function in the continuum limit, their relative overlap has to approach unity. The result 
of such calculation for the interval [0.15 fm, 0.30 fm] is shown in Fig. ||. To understand this, we 
also included two ensembles coarser than £\, with relevant lattice interval being entirely within 
the positive core of the correlator. For S\ the upper limit of the interval involves a correlation 
very close to zero as the correlator just becomes negative. This results in the observed dip and 
large fluctuation of the relative overlap whose definition implicitly assumes that f\ is non-zero on 
its domain (see Eq. (|5])). For finer lattice spacings the relative overlap quickly approaches unity. 
Given these results, we propose to consider the possibility that the following conjecture holds. 

Conjecture 1. Consider the pure-glue lattice gauge theory defined by the Iwasaki gauge action on 
the infinite lattice together with arbitrary non-perturbative procedure of fixing the lattice spacing 
a. IfG(a) and G(a) are the original and the geometric 2-point functions of standard overlap-based 
topological density then 

limO[G(a),G(a), ^,^1 = 1 (6) 
a^o a a 

for arbitrary range of physical distances [r Pt \,r Pi2 \ 

We emphasize that we do not mean to imply that the simplest sign reduction q(x) — ► sgn(g(x)) 
considered here is the unique way to obtain the above-proposed equivalence. Rather, what we wish 
to convey is that there exists a family of computable configuration-based reductions {q(x)} — > 
{sq(x) } with sq(x) G {—1,1} such that the shape equivalence holds. 
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Window: [0.15 fm, 0.30 fm] 
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Figure 2: The relative overlap O[G(a),G(a),0.15/a,0.30/a]. See the discussion in the text for explanation. 

4. The Uses and the Discussion. Dominance of sign-coherent geometry promises to have inter- 
esting ramifications regardless of whether it eventually turns out to be an approximation (exact at 
asymptotically large physical distances) or an exact statement specified by Conjecture 1. In the 
former case, it suggests to consider the "homogeneous approximation" to fundamental topological 
structure [JSj] . In the latter case (strictly consistent with the data at this point) it would represent a 
simple explicit link between short and long-distance physics - a connection that the fundamental 
structure (in any composite field) has to encode in some way. There are several possible ways to 
exploit such exact link, and we mention two examples here. (1) Validity of Conjecture 1 would 
imply that for spectral calculations (glueball masses) it is sufficient to evaluate configurations of 
signs of topological density rather than their precise values. For overlap-based topological density, 
this could lead to a speedup of calculations. (2) An alternative (not strictly equivalent) way of for- 
mulating the observed trends is the following. Consider lim r _,ooG(r,a)/G(r,a). This limit, q 2 {a), 
exists and the correlators G(r,a) and G(r,a) = G(r,a)q 2 (a) are strictly equivalent from the point of 
view of the continuum limit. In particular, in the expected way of obtaining finite physical 2-point 
function we have at arbitrary physical distance r p 

G( r i,a)q 2 (a) 



G P ( r p) 



lim 

a^0 



lim 

a-»0 



lim 

a^0 



(7) 



Applying this equation to short distances x p <C l/m p (glueball range) and taking into account that 
G(r,a) has finite continuum limit at fixed r, one can see that measuring q 2 (a) offers a novel way 
of testing the exact nature of asymptotic freedom. In Fig. || we show the result of evaluating q 2 (a) 
for ensembles S\S^ (fitting G(r,a)/G(r,a) to a constant), indicating that for our finest lattices 
(2.3-3.6 GeV) the behavior is excellently described by q 2 {a) °< a. The perturbative prediction, 
namely a logarithmic decrease of q 2 (a) toward zero, hasn't yet set in. It will be quite interesting 
to examine yet shorter scales to see where the perturbative description starts to dominate, and to 
look for possible geometric changes occurring in the fundamental structure as that happens. The 
absence of such transition would imply a highly non-standard scenario, in which G p {x p ) 1 /\x p \ 7 
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Scale Dependence of qbar 2 




a[fm] 

Figure 3: The scale dependence of q 2 (a).See the discussion in the text for explanation. 

at short distances with the true dimension of the topological density operator being D q = 3.5 rather 
than D q = 4. In that case the parameter 8 would be dimensionful (analogous to mass), asymptotic 
freedom would be only approximate, and the strong CP problem would cease to exist. 1 
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'it is tempting to interpret q(a) as the density of the underlying fundamental structure whose features are fully 
encoded in correlation functions at non-zero distances but distorted by short-range noise at zero physical distance. 
While exact nature of such interpretation depends on the details of the interplay between the noise and the structure, 
it might be sufficiently robust so that the behavior q(a) °= a Da could serve as a definition of the "analytic dimension" 
(£) a -dimensional local density). From Eq. (^|) one could then deduce the relation between the true dimension of the 
operator and the analytic dimension of the associated structure, namely D g +D a = 4. Persistence of the trends shown in 
Fig. H would then imply D a = 1 /2, while the validity of perturbation theory at short distances means that D a = 0. 
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